Melting Scenario for Coulomb-interacting Classical Particles in Two-dimensional Irregular 
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The thermal 'melting' of self-formed rigid structures made of small number of interacting classical particles 
confined in a irregular two-dimensional (2D) space is investigated using Monte Carlo simulations. It is shown 
that the interplay of long-range Coulomb repulsions between these particles and the irregular confinement yields 
a solid-like phase at low temperatures that possesses a bond-orientation order (BOO) but the positional order is 
depleted down to the lowest temperatures. Upon including thermal fluctuations, this solid-like phase smoothly 
crosses over to a liquid-like phase by destroying the BOO. The demise of solidity is aided by the proliferation 
of disclinations. The behavior of different physical observables in the crossover region are obtained. Our results 
will help quantifying melting found in experiments on systems with confined geometries. 



I. INTRODUCTION 

The phenomenon of melting^— has always fascinated 
physicists from the early dawn of the study of condensed 
phases, because it is encountered in everyday life. The or- 
der (or quasi-order) characterizing a solid, as well as the 
mechanism through which a solid melts, depends crucially on 
the spatial dimensionality. 6 For example, a two-dimensional 
(2D) systems present further challenges in comprehending 
melting than in 3D because of the increased significance of 
fluctuations. 7 While serious research efforts have enriched a 
coherent understanding of 2D melting in bulk systems^ the 
physics of melting in confined systems made of small number 
of particles has not been studied as much. While one cannot 
expect a sharp phase-transition 8 in a confined systems as in 
the bulk - even a thermodynamic phase is ill-defined with a 
finite number of particles. However, some signatures of the 
"melting", say, a crossover (CO) from a solid-like to a liquid- 
like phase have been foundi2iI2 

Such 2D-clusters find significant experimental relevance 
and are obtained in a variety of experimental realizations, 
such as, radio-frequency ion traps, 1 1 electrons on the surface 
of liquid He, 12 electrons in quantum dots (QD) in semicon- 
ductor heterostructure, 13 and in dusty plasmas. 14 The inter- 
actions between the charged particles are expected to main- 
tain the bare long range (~ r" 1 ) Coulomb Repulsion (LRCR) 
due to the lack of screening in small clusters compared to 
the bulk. The resulting solid-like phase in a cluster has 
been termed as Wigner Molecule (WM), 13 due to its anal- 
ogy with a Wigner crystal (WC) in bulk system. 15 Inciden- 
tally, a 2DWC undergoes melting to a Fermi liquid even at 
zero Temperature (T = 0), 16 driven by quantum fluctuations 
(QF), and so does a WM in circularly symmetric trapsJ^— 
However, for a variety of experimental clusters at finite T, 
the quantum zero-point motion (QZPM), inherent to particles 
in a bound states, and the associated QF are not so impor- 
tant: either due to the sizable operating T, or because of the 
heavy mass of the particles or their low density, or a com- 
bination of various experimental reasons. The rich physics 
of the CO is then contributed entirely by the classical ther- 
mal fluctuations (TF). T driven melting of classical Coulomb- 
interacting 2D plasma had been studied numerically^ in 
the past, and the corresponding Kosterlitz-Thouless-Halperin- 



Nelson- Young (KTHNY) theory of 2D melting has also been 
developed. 2124 Similar thermal melting of clusters in 2D har- 
monic trap has recently been studied extensively ^ii 2 ! and 
also been realized experimentally, 29 leading to detailed in- 
sight through the measurements of various observables in- 
cluding relative positional fluctuations and associated Lin- 
demann Ratio, 29 pair-correlation function, 30 static structure 
factor 31 and addition spectra. 32 

Progress in the field of high accuracy fabrication made the 
clusters across all different realizations ultra-tunable and the 
shape of confinement is controlled at will by electrostatic and 
magnetic means. 14 These systems naturally have become a 
hotbed for systematic study of the complex interplay of the 
LRCR, quantum interference effects of the confinement, the 
level quantization due to their smallness, and finally the dis- 
order in the form of irregularities on the confinement. Pres- 
ence of irregularities, at least in the large QD have been es- 
tablished in Coulomb blockade experiments, 33 ' 34 and corre- 
sponding theories based on "Universal Hamiltonian"^^ 6 . have 
been put forward. It is this last point - the effect of disorder or 
irregularities on the melting mechanism that has drawn seri- 
ous attention in the recent past 37 in the context of the quantum 
melting 2DWC. 38 While the melting scenario is still some- 
what unsettled even for disordered bulk systems, with propos- 
als of intermediate exotic phases^d . the role of irregularities 
on the melting of WM has received relatively little attention 
(See however, liTl - M) and will be addressed in this work in 
details. 

There are relevant questions with respect to the CO between 
solid-like irregular wigner molecule (IWM) to a liquid-like 
phase: Is there a solid-like phase at all? Does the CO occur 
more or less abruptly or very gradually with Tl What are 
reasonable criteria for identifying the CO in an IWM? And 
finally, what is the physical mechanism behind the 'melting', 
if any? Addressing these important questions within a frame- 
work of systematic calculations, we first show that there in- 
deed exists a 'phase' that is solid-like with the solidity coming 
from the orientational order at low T. Upon inclusion of TF 
induced by T, the IWM crosses over to a liquid-like phase. 
Looking into the behavior of different observables, we esti- 
mate the temperature width ATx for the crossover. We also 
present numerical evidences illuminating the mechanism of 
this melting, which turns out to be associated with formation 
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and proliferation of disclinations on top of the irregular geom- 
etry which already diminished positional order in the IWM. 



II. MODEL AND PARAMETERS 



We consider a system of N classical particles interacting 
via a Coulomb repulsion Vboul = J2i< j (*« 



r 3 y L ( The 

standard LRCR factor C = q 2 /eis assumed unity). The effect 
of irregularity is introduced through the well studied 45 ^ 6 . 2D 
quartic confinement potential V(x, y): 



V c (x,y) 



2Xx 2 y 2 + j(x — y)xyr 



(1) 



where r — ^ x 2 + y 2 . The parameter a makes the confine- 
ment narrow or shallow, we choose b = 7r/4 to break x-y sym- 
metry, A controls the chaoticity, and 7 breaks the reflection 
symmetry. The numerical values of these parameters would 
be scaled by C appropriately to get their dimensional values. 
The classical dynamics of non-interacting particles in above 
confinement have been studied in details 45,46 with the identifi- 
cation of integrable to chaotic regime. The interplay of LRCR 
between spin-1/2 fermions and irregularities of V c (x, y) have 
also been looked atj 4 -^ We choose A € [0.565, 0.635] and 
7 G [0.1,0.2] for our study and generate 5 different realiza- 
tions, each for a given combination of (A, 7) and study the 
behavior of N = 10 — 151 particles and the thermal evolution 
of certain observables. While the particles have no dynam- 
ics at T — due to the lack of QZPM, many low-energy 
configurations with varied position of particles contribute to 
the partition function and hence to the physical observables 
with appropriate Boltzmann weight for T ^ due to the in- 
terplay of thermal kicks on these classical particles and their 
potential energy of LRCR and also from the external irregular 
confinement. The thermal evolution of the system is studied 
using standard Monte Carlo simulated annealing (MCSA) 49 
aided by metropolis algorithm. 50 MCSA is expected to trace 
the appropriate low-energy states consistent with the Boltz- 
mann probability at a given T (T includes ks, the Boltzmann 
factor for all our results), and can obtain the Ground state con- 
figuration in the limit T — >• 0, for an appropriate choice of a 
slow annealing schedule. 51 All physical observables are aver- 
aged over the 5 different realizations for better statistics. 



III. RESULTS 

We begin discussing the results by showing evidences for 
the 'solidity' of the proposed IWM found at T = 0, an exam- 
ple of which is shown in Fig. 1(a) for a fixed set of param- 
eters (See the caption). Such results were obtained by run- 
ning MCSA up to a very low T ~ 10 -5 , and then letting the 
system relax to T = configuration following the downhill 
moves alone in the energy landscape. While it is impossi- 
ble to guarantee accessing the true ground-state configuration 
(GSC) for each individual runs, we ensure that the GSC is in- 
deed achieved by the following statistical analysis: We start 
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FIG. 1. (a) The GSC corresponding to V c (x,y) of Eq.(l) with 
A = 0.2, 7 = 0.635, N = 148 showing the amorphous solid- 
like behavior. Snapshot of particles for the same realization is shown 
in (b) at T — 0.015 where the melting has just started by produc- 
ing a tortuous path connecting the particles in the bulk, and (c) at 
T = 0.065 where the melting is nearly complete. Snapshots of MC 
configurations are presented with intervals in which particles would 
have moved same distance diffusively in both (b) and (c). For com- 
parison, similar snapshot is presented in (d) for hard-core particles 
(core size ~ 0.285) in the same confinement, showing "pre-melting" 
primarily occurring near the boundaries. 



the T = calculations from ~ 200 different initial configu- 
rations generated in the equilibrium MCSA runs at the lowest 
T, for a given set of V c -parameters. Not all these 200 runs 
converge to the same final GSC, however, for all the cases we 
study over the whole parameter space, more than 60% of them 
do with exactly same energy Eq. Furthermore, runs that don't 
converge to the GSC always had E > Eq, yielding us con- 
fidence of identifying the true GSC. We also ensure by using 
first principle Newtonian mechanics that the force, from inter 
particle LRCR and confinement, on each individual particle in 
the GSC adds up to zero, modulo a numerical tolerance. 

Snapshots of MC configurations at T = 0.015 and 0.065 is 
presented in Fig. 1(b) and (c), illustrating the crossover from 
solid-like to liquid-like behavior in the IWM. In Fig. 1(b), 
where the 'melting' has just about started, the thermal mo- 
tion of some particles becomes correlated leading to 'incipi- 
ent' melting through the formation of a tortuous paths. While 
such paths always occur at similar T for all the 5 realizations 
of confinement signaling the commencement of melting, their 
locations are completely random in space without any prefer- 
ence for bulk or boundary. We will bring back our attention 
to these interesting structures when we discuss the possible 
mechanism of the crossover. Fig. 1(c) illustrates the 'melted' 
state, where every particle travels everywhere in the system. 
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Fig. 1(b) should be contrasted with Fig. 1(d) where a similar 
snapshots depicting incipient melting is presented for hard- 
core particles, 52 with average core size same as our Coulomb- 
interacting particles as inferred from pair correlation function 
(defined below), shows that the melting starts primarily on 
the boundaries. While pre-melting on the boundary in solids 
with short-range interactions has been discussed for a long 
timer^ 5 ^ our results with LRCR show that it can occur any- 
where in the systems under considerations. 

Having established that we do track the true GSC, the next 
question is: Are these IWM equivalent to the solid-like phase 
indeed in confined systems? Even the bulk 2DWC leads to 
very interesting physics^^ 5 - Extensive research 56 confirmed 
that an ideal 2D WC phase is characterized by both positional 
and orientational orders, and the demise of both are neces- 
sary to cause its melting, which are likely to occur in two 
stages. 5 ^ Fig. 1(a) clearly shows that the positional order 
(PO) is already depleted even at T — 0, because of the break- 
ing of translational invariance by the confinement leading to 
an amorphous solid-like phase. Lack of positional order is 
also found by calculating p^ — N^ 1 J^i exp(ifc.ri), ft be- 
ing position of i th particle, p^ shows only broad humps in 
our case even at T = and shows no well defined Bragg- 
peaks. This observation points towards absence of PO down 
to T = 0. We also find no significant thermal evolution of p^. 
The orientational order, on the other hand, is evident from the 
nearly perfect 6-coordinated environment for all the particles, 
except obviously for those on the boundary. Our results also 
indicate that a self-formed IWM makes a more or less triangu- 
lar lattice (though with irregular lattice lines in the absence of 
PO) - the typical 2D Bravais lattice minimizing total energy 
with long-range interactions. The quantitative estimate of this 
orientational order is measured through Bond-orientation or- 
der (BOO),— ipQ(r) (Fig. 3), and the Bond-orientated correla- 
tion function (BOCF), g§ (r), to be defined in the next section, 
(Fig. 4), whose evolution with T illustrates the CO from solid- 
like to liquid-like behavior in IWMs, will be discussed in fol- 
lowing sections. The other signature of solidity, even for our 
amorphous solid-like phase, lies in its rigidity leading to the 
absence of any significant positional root mean square (RMS) 
fluctuations of the constituent particles from their equilibrium 
positions at low T. Such fluctuations measured in terms of 
the average inter-particle spacing is known as Lindemann Ra- 
tio (LR), £. 25 Our results in Fig. 2 demonstrates that the av- 
erage LR is essentially zero for all T < 0.01 in the whole 
parameter-space studied, and thus provides further evidence 
for the solid-like IWM at T = 0. 

How does a IWM melt as T is increased? Insight could be 
derived from the LR which we define as: C = N^ 1 J2i £<> 
where d ■ 



-r< 0) . 



i 1 V\K^P) f° r *" tri particle, and u,; = f$ 
Here rf' is the position of the i-th particle in the initial config- 
uration (i.e., the one at the end of the equilibration MC steps), 
and di is the distance of the i-th particle with its neighbors av- 
eraged over all MC configurations. We track the thermal evo- 
lution of £ in Fig. 2. While £(T < 0.01) « 0, LR increases 
dramatically for 0.01 < T < 0.05, beyond which the growth 
of £ becomes far more gradual. This range thus identifies 
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FIG. 2. The main panel shows the evolution of £ (defined in the text) 
with T for different N. While £(T < 0.01) « 0, its smooth rise for 
0.01 < T < 0.05 indicates 'melting'. £(T > 0.06) rises but with 
a much weaker slope, as expected from diffusive motion of particles 
in a 'fluid' phase. The top-left inset presents the T-evolution of the 
probability distribution of d, P(Ci) for TV = 141. The variance 
of P(Ci), ac, (scaled by a factor 155 for clarity) is also presented 
in the main panel, showing similar rise as £ itself. The bottom-right 
inset describes the same crossover as in main panel, but in terms 
of modified Lindemann parameter T (defined in the text), and its 
variance or (scaled by 1200) for TV = 141. 



the crossover width, ATx, between the IWM and its melted 
state. Distribution P(Ci) of £,-, collected over all particles, 
and also over five different realizations of irregular confine- 
ments (shown in the top-left inset of Fig. 2) indicates that a 
peak starts appearing for non-zero values of £; for T > 0.01 
and the variance, ac, of this distribution follow more or less 
similar T-evolution as the £ itself (See main panel of Fig. 2), 
giving further confidence in identifying ATx ■ In the bulk 2D 
systems, £ is difficult to deal with because it can show loga- 
rithmic divergence 26 with system size, and a modified Linde- 
mann parameter, T, is proposed 27 in terms of the relative inter 
particle distance fluctuations: 



r = 



I x N 1 N b 




(2) 



where, the sum on j runs over Ni, number of nearest neighbors 
of z-th particle. While £ cannot show any divergence in IWM 
in contrast to the bulk case due to the finiteness if IWM we, 
nevertheless, calculated the evolution of V and the variance 
or of P(Ti), defined exactly in the same way as those for £^. 
The resulting T and or is shown as a function of T in the 
bottom-right inset of Fig. 2 for N = 141. As expected, £ and 
r show similar behavior including same value of ATx- 
While the thermal evolution of T shows a sharp transition 
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FIG. 3. The distribution, P(ipe), for TV = 125 shows a strong peak 
at i/>6 ~ 1 signifying nearly regular 6-coordinated IWM at low T. 
P(4>6) broadens with T across ATx, finally producing a very broad 
distribution ranging between and 1 beyond the crossover from 
solid-like to liquid-like phase. The boundary-particles, not having 
6 neighbors surrounding them, were not considered in P(ij)§). 

(i.e., AT X -> 0) for the bulk 2D system, 2 ! a finite AT X mea- 
sures the crossover temperature width in a finite system, we 
found that its value for our IWM is typically larger than that 
for a circularly symmetric WM. 28 This can be crudely taken 
as the smearing of ATx by irregularities. Our study also in- 
dicates that ATx does not change much for for N > 35, and 
does not seem to depend on the model parameters [7, A]. Be- 
yond T > 0.05 the thermal evolution of C follows a weaker 
linear trend, as expected for a liquid with particles moving 
around diffusively. 57 T-dependence of LR had been used to 
identify melting in the bulk (having translation symmetry in 
all directions), 27 as well as in circular dots. 17,28 In the latter 
case, the azimuthal symmetry in confining potential and the 
lack of it in the radial direction helps defining radial and az- 
imuthal LR separately, however, we don't have it because our 
confinement breaks all symmetries. 

As mentioned already, GSC for low T in an IWM seems to 
possess 6-coordinated BOO. It is important that we study the 
temperature dependence of the local 6-fold bond-orientation 
order parameter, 58 defined as: 

1 6 

4(r>-^e X p(64) (3) 

fc=i 

Here r is the position of the particle in question, and 9k is the 
angle that the particle at r makes with its six closest neighbors 
k, relative to any arbitrary reference axis. ipe(r) measures lo- 
cal orientational order in the IWM based on the principle that 
all bonds in a perfect triangular lattice should have the same 
6k, modulo 7r/3, implying tpg = 1 for all r. We present the 
distribution P(ipe) m Fig- 3, for N = 125. The sharp peak in 
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FIG. 4. Evolution of BOCF, g 6 (r) for TV = 125, with T, show- 
ing clear Bragg-type peaks at low T signifying bond-oriented order. 
These peaks weaken with increasing T, finally leading to only couple 
of initial humps like a liquid. Data for different T are shifted verti- 
cally for clarity. The boundary particles are not taken into account for 
g&(r), because they don't have 6 surrounding neighbors. Inset shows 
the T-evolution of PVR (defined in the text). The smooth decrease 
of PVR over ATx is indicative of a smooth CO. 



P(ipe) at "06 = 1 demonstrates the persistence of strong BOO 
in the GSC of an IWM at low T, particularly for large N. In 
a liquid, on the other hand, lack of BOO forces Ok € [0, 2tt] 
randomly, resulting into a broad distribution of P(ip6) rang- 
ing between and 1. Our results of T-dependence of P(tl>e) 
indicates that while the change from solid-like to liquid-like 
crossover with T takes place continuously, the most dramatic 
changes occur in same range of ATx where C and ac was 
changing most rapidly (See Fig. 2), furnishing further sup- 
port for identifying ATx as the CO width. For clusters with 
TV < 35, ratio of the number of particles on the boundary 
to those in bulk increases resulting into fewer data-points for 
constructing P(tpQ), augmented with multiple peaks in P(tpe) 
for certain N values even at T = make the distribution qual- 
itatively non-universal. Nevertheless, a 6-coordinated bulk 
particles can be broadly identified down to N ~ 20, at the 
lowest temperatures. The T-dependent broadening of already 
broad P(tpe) for small N does not allow us to identify ATx 
as clearly, and this width is typically larger than that for large 
N. 

While ipe(r) describes strength of local BOO, the long- 
range nature of this order is best discussed in terms of the as- 
sociated bond-oriented correlation function (BOCF), defined 
as: 

g 6 (r) = g 6 (\f\) = (^(1^1)^(1*" - (4) 

which measures the distance up to which local 0g are cor- 
related. It is exactly the same way pair distribution function 
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FIG. 5. Evolution of PDF, g(r) for N = 125, with T, showing much 
weaker Bragg peaks at low T than in ge(r) in Fig. (4), indicating 
weaker positional order than bond-oriented order. The other features 
follow the trend of ga(r). The inset showing a smooth decrease in 
the PVR of g(r), and it is less drastic than that of ge(r) in Fig. (4). 



(PDF), g(r), defined as: 



1 (r)=g 6 (\r\) = {5(\r*\)5(\r*-r\)), 



(5) 



FIG. 6. The main panel shows the distribution P(Ag) for N = 141 
(defined in the text). Traces of individual gi(r) for two particles (1 
and 7 in our nomenclature) are shown on the left-insets for low T, 
showing significant particle-to-particle fluctuations in it, attributed 
by the irregular confinement. At large T, on the other hand, such 
fluctuations vanish, producing gi(r) similar to g(r) for all i, as 
shown on the right-inset for the same two particles as before. The 
resulting P(Ag), thus becomes narrow and symmetric with T, indi- 
cating melting in the same range as ATx as from other criteria. 



measures the positional correlation, describing the probability 
of finding another particle at a distance r away from a given 
one, on the average. Our results for BOCF in Fig. 4 at low 
T shows well defined peaks at specific values of r for up to 
the largest distance between particles in IWM indicating that 
BOO is long-ranged (in the length scale of the IWM). Unlike 
the bulk system, where g§ (r) (and g(r) as well) reaches a con- 
stant for r — >• oo indicating a uniform density, it vanishes in 
IWM beyond the system size because the distance between 
two particles is limited to the linear dimension of IWM. With 
increasing T, the long-range nature diminishes by washing 
out the peaks progressively at large r, and beyond T > 0.05, 
only liquid like behavior persists showing only first couple of 
smeared peaks or humps, as presented in Fig. 4. In the liquid- 
like phase, ge(r) (and g(r) as well) is expected to follow the 
profile of the average radial density for r larger than the po- 
sitions of initial humps, which yields the smooth 7--dependent 
background. The window of T, in which peaks at large r of 
ge{r) starts disappearing, leaving only liquid-like features, is 
consistent with previous ATx- An obvious quantitative mea- 
sure of how fast these Bragg-type peaks are depleting is the 
peak-to-valley ratio (PVR) - the ratio of the value of g§ at its 
first peak at the lowest r to the same on the valley immediately 
after, is also presented as an inset in the same figure. 

We have also studied the thermal evolution of g(r), which 
is presented in Fig. 5. The qualitative features are similar to 
those of g§{r), except that the peaks in g(r) are less sharply 
defined. It is known 5 that for 2D bulk system the long-range 



nature of ge(r) survives up to larger T, than for g(r), lead- 
ing to a hexatic phase 22 at intermediate temperatures. A reli- 
able identification of such a phase from numerics is difficult 
in the IWM owing to its small linear dimension. More in- 
sight on the CO, however, is be obtained from PDF of indi- 
vidual particle, g%{r), in the following manner. gi(r) defines 
PDF of i-th particle with respect to all other (N — 1) parti- 
cles as an average over the MC configurations, in fact, g(r) 
is only the average over gi (r) for all i ranging from 1 ■ ■ • N. 
At T = the unique GSC determines <ft(r) or g(r), unlike 
T 0, where configurations contribute to physical quanti- 
ties with appropriate Boltzman weights, and gi (r) differs a lot 
from site to site due to the irregular confinement depending 
on where in space the i-th particle is located. Thus, the dis- 
tribution, P(Ag), of Ag(r) = gi{r) — g(r) collected for all 
r, and also over different realizations of irregularity, at low T 
will be rather broad, and need not even be symmetric. The 
details of any asymmetry in P(Ag) for a given realization of 
V c (x, y) depends on its parameters, and can be washed out in 
the statistics over many realizations. On the contrary, particles 
are equally likely across the whole system at large T, leading 
to a symmetric and narrow P(Ag). It is then expected that the 
thermal evolution of P(Ag) must shed some light on ATx- 
Note that a similar argument for the bulk system would lead to 
symmetric and possibly narrow P(Ag) for all T. We present 
the T-dependence of P(Ag) in Fig. 6, and the corresponding 
ATx appears to be consistent with other previous criteria. 
Having seen a consistent ATx for solid-like to liquid-like 
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FIG. 7. The specific heat, cv, as calculated from energy fluctuations 
in MC (see text), showing a hump characterizing a crossover in crit- 
ical phenomena. The hump gets sharper as TV increases. Also, it is 
interesting to note that the hump occurs at Tx = 0.03, and is insen- 
sitive to N. The inset compares the TV-dependence of Tx (obtained 
from the location of the hump in cv) and T c , the transition tempera- 
ture for a bulk system with the same average density of particles, as 
in our irregularly confined system for each/V. Modulo the large un- 
certainty in Tx (Note ATx ~ Tx), their behavior indicates that the 
TF could destabilize an IWM more than an equivalent bulk system. 

behavior in an IWM, we turn to the T-dependence specific 
heat (c v ) of the IWM,^ With E = E MC - the total MC 
energy, and E = (E), cy is defined as, 

° V = ^ = {{E2) {E)2)/T2 (6) 

The temperature evolution of cy is presented in Fig. 7, show- 
ing a distinct hump as expected in a crossover, and the hump 
gets sharper as TV increases. It is interesting to note that the 
position of the hump, Tx ~ 0.03, is fairly insensitive to N, 
and falls close to the midway of ATx- The estimate of the Tx 
from the peak-value of smoothed cy is presented in the inset 
of Fig. 7 as a function of N. We also present for compari- 
son the corresponding T c for bulk using VvrW/2c 137 from 
Ref. |20j where for the average density, n, we used its value 
for our IWM from the MC runs. Interestingly, Tx is found to 
lie lower than T c for N > 35. This is qualitatively consistent 
with what happens in a circular Wigner molecule. 28 However, 
serious significance might not be associated with this compar- 
ison, because the uncertainty in Tx is large (ATx ~ Tx)- 
Based on all different criteria for Tx and more broadly for 
ATx that we report here, we find that melting is smeared 
strongly, with the width of transition comparable to the melt- 
ing temperature. 

Our results evidently raise the next important question: 
what's the mechanism for the CO seen here? In the absence of 
a firm analytical answer for the 'melting' in irregular confine- 




FIG. 8. The Voronoi diagram corresponding to single equilibrium 
MC configuration with N = 148 corresponding to Fig. (1). The 
location of the actual particles are shown by thin dots for those with 
CN = 6, filled square for CN = 7, and filled circles for CN = 5. 
Note that CN = 5, 7 implies local disclinations. (a) VD for GSC 
at T = shows predominantly regular hexagonal area surrounding 
each particle, resulting strong BOO. (b) The configuration at T = 
0.015 illustrates the beginning of formation of correlated path of free 
disclinations. (c) The proliferation of disclinations is demonstrated 
on the VD at T = 0.065 (d) The VD for hard-core particles implies 
that the disclinations are limited primarily near the boundary leaving 
sizable BOO in the bulk. 



ments, we turn to numerical evidences and we find the follow- 
ing: As discussed before, melting of 2DWC is mediated by 
production of crystal defects, like dislocations and disclina- 
tions leading to breaking of positional and orientation orders 
respectively. Our results already indicated that the positional 
order is largely depleted in the IWM even at T = and this 
is consistent with our finding of only occasional free disloca- 
tions - which are like an extra row of particles stuck partly, up 
to the lowest T, and more importantly we rarely find any trend 
for T-dependent proliferation of them. Disclinations, charac- 
terized by a mismatch in orientation as one circumnavigates 
it, signal the loss of BOO and are indeed found to proliferate 
significantly for all T > 0.01. They are best seen as a par- 
ticle having the incorrect coordination number as measured 
by the Voronoi diagram (VD) 60 as we discuss below. VD as- 
signs a 2D polygon around each individual particle in a given 
configuration, such that, any point within that polygon will be 
closer to that given particle than others and thus determines 
coordination number (CN) and identifies closest neighbors in 
an unbiased manner. VD corresponding to a perfect triangular 
lattice would be regular hexagons of same size touching each 
other with lattice points lying at the centers. 

We present 3 VDs on an arbitrary equilibrium MC config- 
uration in Fig. 8 for a fixed set of parameters (See the cap- 
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tion), but for three different temperatures. While the results 
at T = shows CN= 6 for nearly all particles (except for 
those on boundary), the VDs at T = 0.05 and 0.1 shows pro- 
liferating disclinations (i.e., CN= 5 or 7) with T, and simi- 
lar qualitative pictures have been found for other parameters 
as well. We now come back to the question: Why does LR 
show sudden rapid increase for T > 0.01? Our numerical 
results indicate that for T > 0.01 (and more or less indepen- 
dent of N) disclinations start growing causing re-adjustments 
of particles in space which expectedly raises C Such re- 
adjustments involve the particles close by, leading naturally 
to correlated movements of particles in certain regions which 
not only explains the incipient melting through the tortuous 
paths found in Fig. lb, but also the large variance ac (Fig. 2) 
for T e ATx as well. Our results indicate that the melting of 
an IWM is associated with disclinations destroying BOO. 

IV. DISCUSSION 

We found that while IWMs with N > 35, present quali- 
tatively similar physics discussed so far - the crossover from 
solid-like to liquid like behavior is lot more smeared, and is 
rife with larger fluctuations in all physical observables stud- 
ied. Further, the ratio of boundary to bulk particles being 
higher, it is harder to look for 'universal' features in those. 
The melting in case of hard-core particles (but no Coulomb 
interactions) commences on the boundary, however, shows 



otherwise similar behaviors for physical observables as our 
Coulomb interacting particles though the scale of ATx is 
larger. 

In conclusion, we have studied the thermal crossover from 
a Wigner-type solid-like to a liquid-like phase in an irreg- 
ular confinement containing classical particles. Our results 
demonstrates that while such crossover takes place gradually 
without any sharp changes, and hence the width ATx is large. 
Thermal evolution of different observables point towards a 
unique ATx- Interestingly, the mechanism for melting ap- 
pears to be the proliferation of disclinations, that destroys 
quasi-long range BOO. Breaking of all symmetries does not, 
as such, stabilize the quasi-order state than bulk. Experiments 
study some of the dynamical properties of the melting, such as 
measuring diffusion constant, frequency dependence of struc- 
ture factor, the dynamics of the defects and their role in melt- 
ing etc. While our method constrains us to study only the 
static properties, an extension to include dynamics of the de- 
fects seems to be a bright direction. We hope that our finding 
will help understanding the physics of chaotic quantum dots 
in the experimental regime. 
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